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■ 

Both the Mie coefficients and the electromagnetic field distributions in the two-layered sphere 
1 containing left-handed media irradiated by a planar electromagnetic wave are presented. 

(SJ ; PACS: 78.35,+c 

x 
S: 

q 1 I. INTRODUCTION 

(N '■ 

More recently, a kind of artificial composite media ( the so-called left-handed media ) having a frequency band where 

■ the effective permittivity ( e ) and the effective permeability ( fi ) are simultaneously negative attracts attention of 
Cl< many researchers in various fields such as materials science, condensed matter physics, optics and electromagnetism 

■ [1-5]. Veselago 1 first considered this peculiar medium and showed that it possesses a negative index of refraction [6]. 

' It follows from the Maxwell's equations that in this medium the Poynting vector and wave vector of electromagnetic 
wave would be antiparallel, i. e., the vector k, the electric field E and the magnetic field H form a left-handed 
system; thus Veselago referred to such materials as "left-handed" , and correspondingly, the ordinary medium in which 

■ k, E and H form a right-handed system may be termed the "right-handed" medium. Other authors call this class 
of materials "negative-index media (NIM)", "double negative media (DNM)" [5] and Veselago's media. It is readily 
verified that in such media having both e and /i negative, there exist a number of peculiar electromagnetic properties, 
for instance, many dramatically different propagation characteristics stem from the sign change of the group velocity, 
including reversals of both the Doppler shift and the Cherenkov radiation, anomalous refraction, modified spontaneous 

1 \ emission rates and even reversals of radiation pressure to radiation tension [2]. In experiments, this artificial negative 

■ electric permittivity media may be obtained by using the array of long metallic wires (ALMWs), which simulates the 
plasma behavior at microwave frequencies, and the artificial negative magnetic permeability media may be built up by 

£SJ using small resonant metallic particles, e. g., the split ring resonators (SRRs), with very high magnetic polarizability 
00 ; [3,7-9]. 

The extinction properties of a sphere ( single-layered ) with negative permittivity and permeability is investigated by 
Ruppin. Since recently Wang and Asher developed a novel method to fabricate nanocomposite SiC>2 spheres ( ~ 100 
nm ) containing homogeneously dispersed Ag quantum dots (2^5 nm ) [10], which may has potential applications 
to design and fabrication of photonic crystals, it is believed that the absorption and transmittance of double-layered 
sphere deserves consideration. In this paper, we present Mie coefficients of double-layered sphere and consider the 
scattering problem, including the topics on field distribution, electromagnetic cross section, extinction spectra as well 
* ^ ■ as some potential peculiar properties arising from the presence of left-handed media. The formulation presented here 
can be easily generalized to cases of multiple-layered spheres. 



X 



II. METHODS 

It is well known that in the absence of electromagnetic sources, the characteristic vectors such as E, B, D, H and 



Hertz vector in the isotropic homogeneous media agree with the same differential equation 



*This paper will be submitted nowhere else for the publication, just uploaded at the e-print archives. 
^E-mail address: jqshen@coer.zju.edu.cn 

x Note that, in the literature, many authors mention the year when Veselago suggested the left-handed media by mistake. 
They claim that Veselago proposed the concept of left-handed media in 1968. On the contrary, the true fact is as follows: 
Veselago's excellent paper was first published in July, 1964 [Usp. Fiz. Nauk 92, 517-526 (1964)]. In 1968, this original paper 
was translated into English by W. H. Furry and published again in the journal of Sov. Phys. Usp. [6]. 



1 



VV-C-Vx VxC + /j 2 C = 0. (2.1) 
The three independent vector solutions to the above equation is [11] 

L = W, M = V x aip, N = ^VxM (2.2) 

with a being a constant vector, where the scalar function tjj satisfies V 2 ip + k 2 ip = 0. It is verified that the vector 
solutions M, N and L possess the following mathematical properties 

M = Lxa=-|-VxN, L • M = 0, VxL = 0, V • L = V 2 tp = -k 2 tp, V • M = 0, V • N = 0. (2.3) 

Set M = mexp(— iwt) and N = nexp(— iut), the vector wave functions M and N can be expressed in terms of the 
following spherical vector wave functions [11] 

m • dP m (cosO) 

me m „ = T^~;Z n (kr)P™ '(cos 0)™m<£i 2 ~ z n (kr) — -™ s m^i 3 , 

sin u ou 

77 f 77 _|_ 1 ) 1/9 f) 

= -4^^(^) P n ro (cos^)^><Aii + ^[rz n (kr)]-P™ (cos ey™mcbi 2 

T l^^WW™ 9 )"™ 4 **- (2 ' 4) 
In what follows we treat the Mie coefficients of double-layered sphere irradiated by a plane wave. 

A. Definitions 

We consider a double-layered sphere with interior radius a\ and external radius a 2 having relative permittivity 
(permeability) t\ (/ii) and e 2 ( fJ- 2 ), respectively, placed in a medium having the relative permittivity eo and perme- 
ability /io- Suppose that the double-layered sphere irradiated by the following plane wave with the electric amplitude 
Eo along the z-direction of Cartesian coordinate system [11] 

2n + 1 



E ; = a x E cxp(ik a z - iwt) = E Q exp(-iwt) ^ i n " (™£in ~ ir 4il) . 

n=i n \ n + i ) 

Hi = a y —E cxp(ik z - iujt) = — —E exp(-iwi) V i n 2 " + ] (n$ n + in^ n ) , (2.5) 



where 



m S" = ±^™( fc or)P^cos0)™ s 0i 2 ~ J„(fcor)^^>i 3 , 

4i« = ^^^(fco^P^cos^^^ir + ^ [fc 0„(fcoO]' ^jf ± ^|n^ [*brjn(*br)]' P^cos0»i 3 (2.6) 

with the primes denoting differentiation with respect to their arguments and k n being ^/eo^o^- In the region r < cl\, 
the wave function is expanded as the following series 

2n+ 1 



E t = S exp(-i W t) ]T ^^yr - <n« ) , 

n— 1 ^ ' 



' n— 1 ^ ' 

Note that here the propagation constant in m^ n and which have been defined in (2.6) is replaced with k\, i.e., 
V £ iMi7- I n the region ai < r < a2, one may expand the electromagnetic wave amplitude as 
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^ n(n+l) 

n=l 

H m = 2 -E exp(-iiut) V t" . 

, rcm + 1) 

n=l v ' 



(2.8) 



(3) (3) 

Note that here in order to obtain m eln and n oln , the spherical Bcssel functions j n in the spherical vector function 
and is replaced by the Hankel functions , namely, the explicit expressions for and n^ n can also 

be obtained from (2.6), so long as we replace the spherical Bessel functions j n in (2.6) with the Hankel functions hn\ 
Apparently, here the propagation constant in and n^ n should be replaced with k 2 , i.e., y/ej/li^- 

In the region r > a 2 the reflected wave is written [11] 

2n + 1 



E r = E Q exp(-^i) ]T i"^^ - i^ngj,) , 



^ n(n + l) 

n+ 1 / 
'^-f n(n + l) [Un ""''" 

n=l x ' 

where the propagation vector in and is replaced by ho, i.e., ^/e ^0 7 



n=l v 
(3) Q ^ „(3) :«, 



(2.9) 



B. Boundary conditions 

At the boundary r = a\ , the boundary condition is given as follows 

ii x E m = ii x E t , ii x H m = ii x H t , 



where the unit vectors of Cartesian coordinate system agree with ii x i 2 = 13, i 2 x 13 = ii, 13 x ii = h- 
It follows from ii x E m = ii x E t that 

<%j n (N 2Pl ) + a*hW(N 2Pl ) = alj n (N lPl ), 
[N 2Pl j n (N 2Pl )}' + Nrf™ [N 2Pl h^(N 2Pl )\ ' = N 2 bl [N lPl j n (N lPl )}' , 

where P i = fc oi, Ni = 

In the similar manner, it follows from ii x H m = ii x H t that 

Ml (<C [N 2PlJn (N 2Pl )}' + < \N 2Pl h^(N 2Pl )]'\ = P2 {a* [JV lPl j n (JV lPl )]'} , 



V%MN 2Pl ) + b™hW(N 2Pl )\ = N lP2 blj n (N lPl ). 

At the boundary r = (Z2, in the same fashion, it follows from ii x (E; + E r ) = ii x E m that 

j n ( P2 ) + a^\ P2 ) = a™j n (N 2P2 ) + a™h£\N 2P2 ), 
N 2 [ P2jn ( P2 )}' + N 2 V n [ P2 h^( P2 )]' = C [N 2P2 j n (N 2P2 )]' + C [N 2P2 hW(N 2P2 ) 

where p2 = k a 2 , N 2 = f^. 

In the meanwhile, it follows from ii x (H; + H r ) = ii x H m that 

M2 {< [p2h£\ P2 )] ' + [ P2jn ( P 2)} j = Mo {<C [N 2P2 j n (N 2P2 )}' + a* [N 2P2 h£\N 2P2 ) 
^ n h^\ P2 ) + ViiM = N 2Po b™j n (N 2P2 ) + N 2Po b^(N 2P2 ). 



(2.10) 



(2.11) 



(2.12) 



(2.13) 



(2.14) 
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C. Calculation of Mie coefficients 



1. Mie coefficient a r n 



According to the above boundary conditions, one can arrive at the following matrix equation 

-3n(N lPl ) Jn(N 2Pl ) h£\N 2Pl ) \ 

-M2 [JVlPljn(JVipi)]' fJ-i[N 2 pijn(N 2 pi)}' Mi [N 2 pih ( n } (N 2Pl ) 
3n{N 2P2 ) h^(N 2P2 ) 








V 







PO [^2P2in(^2P2)]' N 2 p2^ ) (N 2 p 2 ) 



~P2 



(/*) 
P2kn\ P2 ) 



J 



For convenience, the 4x4 matrix on the left handed side of (2.15) is denoted by 

\ 



A = 



f A n A 12 A 13 

A 21 A 22 A 23 

A 32 A 33 A 34 

\ A 42 A 43 A 44 J 



whose determinant is 
det.4 



= A XX A 22 A 33 A 44 - A n A 22 A 34 A 43 - A n A 32 A 23 A 44 + A n A 42 A 23 A 34 
- A 21 A 12 A 33 A 44 + A 21 A 12 A 34 A 43 + A 21 A 32 A 13 A 44 - A 21 A 42 A 13 A 34 . 



So, the Mie coefficient a r n is of the form 



with 



A r = 







detA r 








detyl 




[A n 


A 12 


A 13 





A 2 i 


A 22 


A 23 








A 32 


A 33 


3n{P2) 



\ A 42 A 43 fl 2 [p 2 j n (p 2 )]' J 



whose determinant is 



\ 



/ 


3n(p2) 
V M2 [P23n(P2)]' J 



(2.15) 



detA r = A n A 22 A 33 p 2 [p2j n {P2)}' - A n A 22 j n (p 2 )A 43 - A n A 32 A 23 p 2 [p 2 j n (p 2 )]' + A n A 42 A 23 j n (p 2 ) 
- A 21 A 12 A 33 p 2 [p 2 j n (p 2 )]' + A 21 A 12 j n (p 2 )A 43 + A 21 A 32 A 13 p 2 [p 2 j n (p 2 )}' - A 21 A 42 A 13 j n (p 2 ). 



(2.16) 



(2.17) 



(2.18) 



(2.19) 



(2.20) 



2. Mie coefficient ¥ n 

According to the above boundary conditions, one can arrive at the following matrix equation 



( -N 2 [N lPl j n {N lPl )\ Nx[N iPx j n {N iPx yl Nx\N iPx h&\N2p x ) 
-N lP2 j n (N lPl ) N 2 mj n {N 2 pi) 



N 2 ^h { r l\N 2Pl ) 



V 



[N 2 p 2 j n (N 2 p 2 )}' 
N 2 p jn(N 2 p 2 ) 



N 2 p 2 h£\N 2 p 2 ) 
N 2 p h { n\N 2 p 2 ] 



-No 






P2h ( h\ P 2) 

,(1) 



\ 



-p 2 h\ '(p 2 ) J 



f bt A 

K 
K 



( 



\ 



N 2 [p 2 jn(p2)}' 
V M2.?n(p2) J 

(2.21) 



For convenience, the 4x4 matrix on the left handed side of (2.21) is denoted by 
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B = 



I B\\ B\2 B13 \ 

B21 B22 B23 

B32 B33 B34 

\ B42 -B43 B44 ) 



whose determinant is 



deti? = B11B22B33B44 — B11B22B34B43 — B11B32B23B44 + B11B42B23B34 
— B21B12B33B44 + B21B12B34B43 + B21B32B13B44 — B21B42B13B34. 



So, the Mie coefficient ¥ n is of the form 



bi = 



dct5 r 
dctB 



with 



B r 



/fill 


Bu 


B\3 







B21 


B22 


B23 










B32 


B33 


^2 [P2jn(p2)}' 




V 


B42 


B43 




) 



(2.22) 



(2.23) 



(2.24) 



(2.25) 



whose determinant is 

det-B r = B n B22B33^ 2 jn(P2) - BnB 2 2N 2 [p2.jn{P2)}' B43 - B n B32B23 P2 jn {p2 ) + BHB42B23N2 fajnfa)]' 

- B2\B 1 2B33H23n{p2) + B2lB 12 N 2 [p2jn(p2)}' B43 + B 2 lB3 2 B 1 3fJ,2jn(p2) - B21B42B13N2 [p2jn(p2)]' ■ (2.26) 



III. DISCUSSION OF SOME TYPICAL CASES WITH LEFT-HANDED MEDIA INVOLVED 



In this section we briefly discuss several cases with left-handed media involved. 



(i) If the following conditions hn\N2Pi) = 0, 



h { n\N 2 pi) 



0, h { n 1) (N 2 p2) 



h [ n\p2), h> n L, {N 2 p2) 



h£\ P 2) 



N 2 = P2 = £2 



— 1 and yu = +1 are satisfied, then it is easily verified that 
det^4 = 0, det£> = and ai = oo, UL 



oo. 



(3-1) 



Thus in this case the extinction cross section of the two-layered sphere containing left-handed media is rather large. 



0. 



(ii) If the following conditions hn\N2p\ 
N2 = P2 = £2 = — 1 and fi = +1 arc satisfied, then it is easily verified that 

detA r = 0, detB r = and a r n = 0, b T n = 



^(iV2Pl) =0, h^(N 2 p 2 ) =jn(P2), h£ } (N 2 p 2 ) =\j n (p2)}', 



(3.2) 



Thus in this case the extinction cross section of the two-layered sphere containing left-handed media is negligibly 
small. 

(iii) The case with N\ = fii = e\ = —1, N 2 = p.2 = £2 = 0, N a = /i = eo = +1 is of physical interest, which 
deserves consideration by using the Mie coefficients presented above. 

Based on the calculation of Mie coefficients, one can treat the scattering problem of double- layered sphere irradiated 
by a plane wave. The scattering and absorption properties of double-layered sphere containing left-handed media can 
thus be discussed in detail, which are now under consideration and will be submitted elsewhere for the publication. 



Acknowledgements This project is supported by the National Natural Science Foundation of China under the 
project No. 90101024. 
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Appendix A: Reduced to the case of single-layered sphere 

In order to sec whether the above Mie coefficients in the case of double-layered sphere is correct or not, we consider 
the reduction problem of double-layered case to the single-layered one when the following reduction conditions are 
satisfied: A n = -A i2 , A 2 i = -A 22 Ni = N 2 , p\ = p 2 , pi = fi 2 . 

By lengthy calculation, we obtain 



and 



Thus 



detA = {A n A 23 - A 21 A 13 ) {A 42 A 3A - A 32 A 4A ) 

= (AnA 23 - A 21 A 13 ) ^j n (N 2P2 )n 2 p 2 h£\p 2 ) - p [N 2 p 2 j n (N 2 p 2 )}' ' h£\p 2 ) 

detA r = (A n A 23 - A 21 A 13 ) {p [N 2 p 2 j n (N 2 p 2 )}' j n (p 2 ) - j n (N 2 p 2 )p 2 [p 2 j n (p2)}'} ■ 

r _ MO [N 2 p 2 j n (N 2 p 2 )]' j n {p 2 ) - jn(N 2 p 2 )p 2 [p 2 jn{p2)}' 



jn(N 2 p 2 )p 2 p2h { n \p 2 ) - Pa[N2P2jn(N 2 P2)]' h ( n'(p2 



(Al) 



(A2) 



(A3) 



which is just the Mie coefficient a T n of single-layered sphere [11]. 

In the same fashion, when the reduction conditions Bu = — -B12, -B21 = —B22 Ni = N 2 , pi = P2, Pi = P2 arc 
satisfied, one can arrive at 



and 



So, 



det£? = (£?iil?23 — -B21B13) (B42-B34 — -B32B44) 

= (B11-B23 — B 2 \Bi 3 ) |[A r 2P2j„(A r 2P2)]V2/il 1) (P2) ~ Nlp a j n {N 2 p2) [pik™ (fti) 

dctB r = (B n B 23 - B 2 lB 13 ) {N^p j n (N 2P2 ) [p 2 Jn(p2)}' ~ [N2P2jn(N 2 p2)} ' P2jn{p2)} ■ 
_ Nlpojn{N 2 P2) [p23n{P2)]' ~ [N 2 p23n{N 2 P2)\' P23n{Pl) 



[N 2P2 ] n {N 2P2 )\ ' p 2 h { n ] (p 2 ) -Nlfi j n (N 2 p 2 ) \p 2 ti n L) \p 2 ) 

which is just the Mie coefficient ¥ n of single-layered sphere [11]. 

This, therefore, means that the Mie coefficients of double-layered sphere presented here is right. 



(A4) 



(A5) 



(A6) 



Appendix B: All the Mie coefficients 

In this appendix, all the Mie coefficients in the double-layered sphere are given. The results are as follows: 



The 4x4 matrix A t is 



A. Mie coefficients a l n , a™ and a™ 



/ A 12 A 13 ^ 

A 22 A 23 

3n{p2) A 32 A 33 A34 

V P2 [p2jn{P2)}' A 42 A 43 A 44 J 



(Bl) 



G 



whose determinant is 

dctA t = jn{p2)A 12 A 23 A 44 - j n (p 2 ) A 22 A 13 A 44 - M2 \Pl jn (P2 )]' A 12 A 23 A 34 + M2 [P2j'n(P2)]' A 22 A 13 A 34 . 

So, 

t detA t 
°" = dctX' 

The 4x4 matrix A m is 



/An 


A 13 







A 2 i 


A 23 










jn(p2) A 33 


A 34 




V o 


Pi [p2jn(P2)]' A 43 


A 44 


) 



whose determinant is 

detA m = -Anj n (p 2 )A 23 A 44 + A n fi 2 [p2j n {P2)]' A 23 A 34 + A 21 j n (p 2 )A 13 A 44 - A 21 p 2 [p 2 j n (P2)]' A 13 A 3i . 

So, 

m _ detA m 
°" ~ detA " 

The 4x4 matrix A & is 

/An A 12 \ 

_ A 2 i A 22 

rn " A 32 Jn ( P2 ) i A 34 

\ A 42 p 2 [p 2 jn(P2)}' A 44 J 

whose determinant is 

dctA ffl = A 11 A 22 j n (p 2 )A 44 - A n A 22 A 34 p 2 [p 2 jn(p2)}' - A 21 A 12 j n (p 2 )A 44 + A 21 A 12 A 34 p 2 [p 2 jn(p2)}' ■ 

So, 

& _ detA & 
°" ~ detA • 



It is readily verified that under the reduction conditions An = —A\ 2 , A 2 \ = —A 22 , N\ = N 2 , pi = p 2 , 
one can arrive at 

dctA t = detA m , a\ = a™, detA & = 0, a™ = 0, 
which means that the above Mie coefficient can be reduced to the those of single-layered sphere. 

B. Mie coefficients b n , 6™ and 6™ 

The 4x4 matrix B t is 

/ B 12 B 13 \ 

q B 22 B 23 

* N 2 [P2jn(p2)]' B 32 B 33 B 34 

V P2jn{P2) B 42 B 43 B 44 ) 



whose determinant is 



detBt = N 2 [P2jn(p2)]' B 12 B 23 B 44 - N 2 [p2jn{P2)}' B 2 2B 13 B 44 - p 2 jn{P2)B 12 B 23 B 34 + p 2 jn(p 2 )B 22 B 13 B 34 . (B12) 



So, 



t detgt 
" dctS 



(B13) 



The 4x4 matrix B m is 



B n 



/fill 





B\ 3 





\ 


■B21 





B2 3 










^2 [P23n{p2)\ 


B 33 


-B34 




V 


P-2jn{P2) 


B 43 


-B44 


/ 



(B14) 



whose determinant is 

detB m = -B n N 2 [p 2 j n {p2)}' B 23 B 4A + B 11 p2jn{p2)B2 3 B 34 + B 21 N 2 [p 2 jn{p2)]' -B13B44 - B 21 ^, 2 j n (p 2 )Bi 3 B 34 . (B15) 
So, 

detB m 



dctB 



(B16) 



The 4x4 matrix B rn is 



B 12 ^ 

^21 -B22 

B 32 N 2 [p 2 j n (p 2 )]' B 34 

V B 42 p, 2 jn{P2) B 44 j 



(B17) 



whose determinant is 

detBa = B X1 B 22 N 2 [p 2 j n (p 2 )}' B 44 - B 11 B 22 B 34 p 2 j n (p 2 ) - B 21 B 12 N 2 [p 2 j n (p 2 )}' B 44 + B 21 B 12 B 34 p 2 j n {p 2 ). (B18) 
So, 

det Rt, 



dctB 



(B19) 



It is readily verified that under the reduction conditions Bu = — B i2 , B 2 \ = — B 2 2, N\ = N 2 , pi = p 2 , pi = p 2 , 
one can arrive at 



detBt = det£ m , b\ = 6™, det^ = 0, 6* = 0, 
which means that the above Mie coefficient can be reduced to the those of single-layered sphere. 



(B20) 
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